Resonant cavity photon creation via the dynamical Casimir effect 
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Motivated by a recent proposal for an experimental verification of the dynamical Casimir effect, 
the macroscopic electromagnetic field within a perfect cavity containing a thin slab with a time- 
dependent dielectric permittivity is quantized in terms of the dual potentials. For the resonance 
case, the number of photons created out of the vacuum due to the dynamical Casimir effect is 
calculated for both polarizations (TE and TM). 
PACS: 42.50.Lc, 03.70.+k, 42.50.Dv, 42.60.Da. 
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One of the most impressive manifestations of the non- 
trivial structure of the vacuum is the static Casimir effect, 
i.e., the attraction of two perfectly conducting plates, for 
example, generated by the corresponding distortion of 
the electromagnetic vacuum state [1]. The non-incrtial 
motion of a mirror can even create particles (i.e., pho- 
tons) out of the vacuum [2] due to the time-dependent 
disturbance - which is called (in analogy) dynamical 
Casimir effect (see, e.g., [3] for review). Unfortunately, 
in contrast to the former (static) effect, the latter (dy- 
namical) prediction has not been experimentally verified 
yet. To this end, it is probably advantageous to exploit 
the drastic enhancement of the number of created pho- 
tons within a cavity occurring if the frequency of the wall 
vibration is in resonance with one of the (discrete) cavity 
modes. The difficulty of accomplishing mechanical vibra- 
tions of the wall with high frequencies (and appropriate 
amplitudes) experimentally has lead to the idea of sim- 
ulating the wall motion by manipulating the dielectric 
permittivity (or magnetic permeability) of some medium 
within the cavity (which can be done much faster). E.g., 
filling the whole cavity with a homogeneous medium de- 
scribed by a time-dependent permittivity s{t) is analo- 
gous to introducing an effective length of the cavity via 
Losit) = \/e{t) L. However, since it is rather difficult 
to influence a medium filling the complete cavity, a new 
proposal [4] (see also [5]) for an experimental verification 
of the dynamical Casimir effect envisions a small slab 
with a fixed thickness a and a time-dependent permittiv- 
ity e{t) located at one of the walls of the cavity, cf. Fig. 1. 
The question of whether and how the motion of the cav- 
ity wall can be simulated by such a small dielectric slab 
- especially in view of the number of created photons - 
will be the subject of the subsequent considerations (see 
also [6] for a 1-1-1 dimensional scalar field model). 

In this Letter, we present an ab initio derivation of the 
dynamical Casimir effect based on the quantization of the 
full macroscopic electromagnetic field within a (perfect) 
cavity with space-time dependent dielectric properties - 
superseding previous effectively 1-1-1 dimensional calcu- 
lations (scalar field model, see, e.g., [6,7]) and approaches 
based on special factorization assumptions (see, e.g., [8]). 
For 3-1-1 dimensional cavities with moving walls, there 
exist various calculations for scalar fields, but very few 



taking into account the full electromagnetic field. E.g., 
in [9] , the electromagnetic field is effectively split up into 
two independent scalar fields obeying different boundary 
conditions via introducing different potentials for the TE 
and the TM modes - which leads to a decoupling of the 
polarizations (TE and TM) per construction. However, 
in the most general situation, TE and TM modes can 
mix - hence their coupling should not be excluded a pri- 
ori but investigated for each special case. 
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FIG. 1. Sketch of the (lossless) cavity containing a thin 
slab with a time-dependent dielectric permittivity (t). 

Since we are considering low-frequency (e.g., mi- 
crowave) photons only, we start from the macroscopic 
source- free Maxwell equations {eq = /io = fi = 1) 

VB = VD = 0,D = VxH,B= VxE, (1) 

with H(t,r) = B{t,r) and D{t,r) = e{t,r)E{t,r). If 
we were to use the usual vector potential A in temporal 
gauge ($ = 0) 



A, B = V X A. 



(2) 



the constraint V • [e{t,r)A{t,r)] — would render the 
usual canonical quantization 



A,it,r),Dj{t,r') =Cy(r,r') 



(3) 



in connection with eliminating the longitudinal degree of 
freedom rather tedious (cL also [10]) because, in this case, 
V • £) = implies djCij = but diCij 7^ in general. 

Therefore, we avoid these difficulty with the well- 
known trick of introducing the dual vector potential (see, 
e.g., [9]) 



1 



H = A, D = V X A, 



(4) 



which applies in this form to the source-free Maxwell 
equations (1) only. In terms of the dual vector poten- 
tial, the constraint simply reads V • A = 0. After the 
duality transformation [11], the Lagrangian is still the 
usual Larmor invariant - but with the opposite sign 



£ = ^ j d^r[B-H-E-D] , 
and the Hamiltonian is again the total energy 



(5) 



= 11 



B H + E D] 

A' + i(V X A)2 



(6) 



The continuity conditions (AA = A^ — A^^) for the dual 

vector potential at the interface between the regions / 
and II of the cavity can be derived from the Maxwell 
equations (1) 



AA = A(Vx A)j^ = A( ^Vx A) =0. 



(7) 



Assuming that the walls of the cavity are perfectly con- 
ducting, for example, the boundary conditions read 



(V X A)|| = (A X [V X A])^ = . (8) 



Consequently, the boundary term arising from the inte- 
gration by parts (as in the Poynting theorem) of the term 
(V X A)^ in Eq. (6) vanishes. Hence we can introduce a 
non-negative and self-adjoint operator /C via 



^/a = V X - V X / J = 



(9) 



with eigenfunctions and eigenvalues O^. Note that we 
consider a lossless (ideal) dielectric medium resulting in 
a real permittivity e G M (and hence a self-adjoint oper- 
ator /C). Owing to the time-dependence of the dielectric 
permittivity e{t,r), the operator IC{t) and consequently 
its eigenfunctions fait^f) as well as eigenvalues fi^(t) 
are also explicitely time-dependent in general. 

The longitudinal modes f\ form the (orthogonal) 
eigenspace with zero eigenvalue V x /|[ =0 and hence we 
can restrict the operator /C to the constraint sub-space 
V • = 0. Since /C is a real operator, we can choose its 
eigenfunctions to be real as well = /* ; and because 
JC is self-adjoint, its eigenfunctions are orthonormal (for 
equal times) 



/ 



(10) 



and complete 



(11) 



with d^f (r — r') denoting the transversal Dirac S- 
distribution diS^ir — r') — 0. Hence a corresponding 
normal mode expansion of the Lagrangian and the Hamil- 
tonian in terms of the dual potentials into the instanta- 
neous basis 



A(t,r)=^g„(t)/„(t,r), 



(12) 



leads to (see also [12]) 

From now on, we shall drop the summation signs for 
convenience by declaring a corresponding (Einstein-like) 
sum convention. The canonical conjugated momenta are 
given by Pa = Qa + Mai3{t)Q(3 and the anti-symmetric 
inter-mode coupling matrix reads 



/ 



(14) 



The usual equal-time canonical commutation relations, 
e.g., [QcPff] = i5ai3, are equivalent to the commutators 
for the fields, such as [AJ(t,r),n'=(t, r')] = iS^^ir - r'). 

It will be convenient to classify the cigcnmodes with 
respect to their polarization at the interface into TE 
(transversal electric) and TM (transversal magnetic) 
modes 



E' 



TE 



0, Bi 



TM 



0. 



(15) 



In terms of the dual potentials, these conditions read 
(V X A)'^^ — and A™ — 0. Assuming the absence of 
any static fields and fields outside the cavity (the walls 
are supposed to be perfectly reflecting), the boundary 
condition B± = implies A_l = and hence Eq. (8) 
imposes the condition V_lA|| =0. As a result, we can 
make the following separation ansatz in the homogeneous 
region / 



fk 



h(T 



sin(fc^a;) cos{kyy) cos{k^z)ei^^ 
cos{klx) sm{kyy) cos{k:,z)e{^y 
cos{kix) cos{kyy) sin(fc^^;)e^^_. 



(16) 



and analogously for region II with x being replaced by 

X — L. 

The wave-numbers ky and kz are simply determined 
by the perpendicular cavity dimensions Ly, via ky = 
Uyw/Ly and kz = UzTt/Lz, respectively, with integers 
Uy and Uz- The remaining polarization factors el/J^ as 

well as the fcx^'^^-values have to be determined according 
to the continuity conditions in Eq. (7), the polarization 
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condition (TE or TM) in Eq. (15), the transvcrsality con- 
dition V • A = 0, the overall normalization in Eq. (10), 
and, finally, the eigenvalue equation (for a fixed time) 



^2 _ (^a)^ + fey + 



[ki^f + kl+kl 



(17) 



which provides a relation between A;^ and k^.^ . Using the 
conditions mentioned above, we arrive at the transcen- 
dental equations 



TE : 



TM : 



tan(aA;^) tan(fc^^[a — L]) 
hi ~ hii ' 

fc^tan(aA;^) _ A;^^ tan(A;^^[a - L]) 



(18) 



which have to be satisfied simultaneously to the eigen- 
value equation (17). 

Assuming the slab to be sufficiently small a <C i, we 
can find approximate solutions for the TE modes 



Ji 



IT 



and for the TM modes (for > 0) 



k^ — 



a 
L 



-II 



- 1 



+ 



(19) 



(20) 



with fcjj = fc^ + k1 and k^ = UxT^/L. 

We observe that the first-order (in a/L 1) contribu- 
tions to the eigenvalues of the TE modes are indepen- 
dent of £^/^^{t). Only for the TM modes, a variation of 
the permittivities e^/^^ [t) induces a change of the eigen- 
values (with the label a = {n, TM}) 



Aft 



-II 



1 



O 



J,]- (21) 



The first-order term of the coupling matrix can be de- 
rived in complete analogy: Mafj{t) oc [a/ L)dt{e^^ / e^). 

In order to simulate an oscillation of the wall, we as- 
sume a harmonic time-dependence of the ratio 



-II 



(22) 



with the amplitude x a-nd an irrelevant additive constant 
^ (which just induces a constant shift of the eigcnfre- 
qucncics). A small harmonic perturbation over a rela- 
tively long time duration (i.e., many oscillations) enables 
us to employ the rotating wave approximation, which ne- 
glects all non-rcsonant terms. According to Eq. (13) with 
nl(t) = (ftO )2 -)- Anl{t), the perturbation Hamiltonian 
can be split up into two parts, the diagonal (so-called 
squeezing) term AQ'^{t)Q'^/2 and the off-diagonal (so- 
called velocity) contribution PaQfiMafj{t), cf. [12]. The 
resonance condition for the former (squeezing) term reads 



w = 2f7^ and for the latter inter-mode coupling (veloc- 
ity) contribution w = |r2° ± f2^|. In the following, we 
shall assume a cavity with well-separated eigenfrequen- 
cies where the external oscillation frequency uj matches 
the diagonal resonance condition uj = 2f2° for a certain 
TM mode only (no resonant inter-mode coupling). In 
this case, the effective Hamiltonian reads 



^•^^ ~ 2 a; e^^ 



tl2 



o 



(23) 



for the resonant TM mode a. Accordingly, the time- 
dependence of the dielectric permittivity of the thin slab 
induces the creation of an exponentially increasing num- 
ber of particles (photons) out of the vacuum (dynamical 
Casimir effect) 



{N){t) = (0| eM+i^est} eM-mst] |0) 



(24) 



Let us state the physical assumptions entering the above 
derivation. Firstly, by starting from the source-free 
macroscopic Maxwell equations with perfectly conduct- 
ing boundary conditions we assumed an ideal cavity and 
omitted losses and decoherence etc. Of course, the ap- 
plicability of this assumption has to be checked (e.g., 
whether the Q-factor of the cavity is large enough) be- 
fore conducting a corresponding experiment. Secondly, 
the external oscillation was assumed to be harmonic with 
the frequency matching the resonance condition exactly. 
Other periodic time-dependences could lead to the con- 
tribution of higher harmonics (2u;, etc.) and one has to 
make sure that the possibly resulting inter-mode coupling 
docs not spoil the main contribution in Eq. (24) . A devi- 
ation (detuning) from the exact resonances; — 2r2j^(l+5) 
is also not critical as long as the relative detuning 5 is 
smaller than the relative perturbation amplitude, cf. [13]. 
Thirdly, the neglect of the higher-order terms in the Tay- 
lor expansion of the transcendental matching equations 
(18) leading to Eqs. (19) and (20) assuming a small slab 
o <C L is only justified if all other involved quantities are 
not too large. If the ratio e^^ /e^{t) changes drastically, 
this approximation breaks down as soon as the small- 
ness of the expansion parameter a/L is compensated by 
a huge variation in e^^ /e^{t). 

Let us study the two fimiting cases: For » e^^, the 
wave-numbers behave as k^ ^ k^^ according to Eq. (17). 
Hence the poles of the functions tan(afc^) in Eq. (18) 
induce drastic changes of k^^ and thus ft - for both, 
TE and TM modes. In this case, the eigenmodes are 
'pulled' into the small slab (which is not desirable since 
the time-varying material properties will entail the dan- 
ger of dissipation and decoherence). In the opposite case 
£^ ^ e^^, the wave- number in the small slab becomes 
imaginary, cf. Eq. (17). In some sense, the modes are 
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'pushed' out of the small slab similar to the phenomenon 
of total reflection (as in an optical fiber, for example) in 
this situation. With exactly the same argument as be- 
fore, there is no effect to lowest order in a/L for the TE 
modes. However, the values kl.^ and hence Vl of the TM 
modes change strongly owing to the occurrence of the 
term e" /e\ cf. Eqs. (20) and (21). 

In summary, a (small and smooth) motion of the 
cavity wall (more precisely, its resonant features) can 
only be simulated for TM modes with i™ = L[l — 
{a/L){e"/e^ - l}(fcf|/fci)] + 0{a?/L'^) via moderate 

changes of e^^ /e^{t) ~ huge variations are inappropriate. 
At a first glance, this result might be a bit surprising since 
the initial (e.g., s^^ = e^) and final situations (e^^ and 
vastly different) are similar to two cavities with slightly 
different dimensions. However, one has to bear in mind 
that the whole evolution - and not just the initial and fi- 
nal state - is important for the dynamical Casimir effect. 
E.g., instead of moving one wall of the cavity, one could 
mediate between the same initial and final states by in- 
serting an additional wall into the cavity - which would 
imply a completely different dynamics (e.g., depending 
on what happens in the cut-off part of the cavity). 

With the aid of the dual vector potential A, we were 
able to quantize the (macroscopic) electromagnetic field 
within a cavity with space-time dependent dielectric 
properties e{t,r) and /z = const facilitating the inves- 
tigation of the influence of the polarization (TE and TM 
modes) etc. In the opposite case ^{t, r) and e = const, 
an analogous derivation can be accomplished using the 
ordinary vector potential A. In contrast to the former 
case, where only the TM modes feel (to first order) the 
presence of the dielectric slab, both polarizations (TE 
and TM) arc effected by the changing magnetic prop- 
erties /u(i,r) in the latter situation [14]. The physical 
difference between the two cases can be explained by the 
distinct boundary conditions which involve real macro- 
scopic charges and currents and are therefore not invari- 
ant under the duality transformation [11]. 

In conclusion, there are three major differences be- 
tween the scenario under investigation and a cavity with 
a moving wall: firstly, the dependence on the polariza- 
tion (only TM modes), secondly, the dependence on the 
perpendicular wave- number oc fc|, and thirdly, that fact 
that the effect does not vanish for = - in this case, 
one obtains just half the value given in Eq. (21). 

Let us insert some explicite numbers: If we assume 
a switching time of about 100 ps [4], the term in 
Eq. (24) is of order GHz. According to the resonance 
condition, the dimensions of the cavity and hence the 
wave-length of the created photons should be several cen- 
timeters (i.e., microwaves). The second term x/s^^ can- 
not exceed 1/2 and is assumed to be of order one. The 
remaining parameter a/L should be small, say 1/100 - 
where the explicite value of a depends on the way of 



changing the dielectric properties, e.g., laser illumination 
depth [4]. Note that a/L = 1/100 is stiU by far larger 
than the relative amplitudes that can be achieved by me- 
chanical vibrations (at these frequencies) . Provided that 
the assumption of a perfect cavity (e.g., Q-factor) is ap- 
propriate during such time-scales, one would create a sig- 
nificant amount of photons after a few microseconds. 
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